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Abstract. We present a 4.5 — 4.85 /im R — 5 000 spectrum of the low mass class I young stellar object GSS 30 
IRSl {L = 25 Lq) in the p Ophiuchus core, observed with the infrared spectrometer (ISAAC) on the Very Large 
Telescope (VLT-UTl). Strong line emission from the ro-vibrational transitions of '^^CO and ^"^CO is detected. 
In total more than 40 distinct lines are seen in the covered region. The line emission is spatially extended and 
detected up to 2" = 320 AU from the central source but is spectrally unresolved (Aw < 30 km s~^). This is the 
first time strong emission in the fundamental ro-vibrational band from CO has been observed from an embedded 
young stellar object. The line fluxes were modeled using a 1-dimensional full radiative transfer code, which shows 
that the emission is fully consistent with a gas in LTE at a single well constrained temperature (T — 515 ± 5 K). 
Furthermore, the ratios between lines from the two detected isotopic species of CO show that the ^^CO lines 
must be optically thick. However, this is inconsistent with the observed spatial extent of the emission, since this 
implies such low CO column densities that the lines are optically thin. A likely solution to the discrepancy is that 
the lines are emitted by a smaller more dense region and then scattered in the bipolar cavity present around the 
central star. This gives a rough estimate of the total molecular gas mass of 1 — 100 and a physical extent 
of ~ 20 — 100 AU. We propose that the most likely origin of the line emission is post-shocked gas in a dense 
dissociative accretion shock from the inner 10 — 50 AU of a circumstellar disk. The presence of a shock capable 
of dissociating molecules in the disk will have implications for the chemical evolution of disks around young low 
mass stars. 
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1. Introduction 

The innermost 50 AU of the circumstellar environ- 
ments around embedded young stellar objects is poorly 
constrained observationally due to the large extinction 
through the embedding material {Ay > 20 mag) and the 
small angular size (< 1") of the region. The usual molec- 
ular probes in the millimeter-submillimeter region are not 
effective for the temperatures (T > 200 K) and densities 
(riHa > 10'' cm~^) thought to be present. However, an un- 
derstanding of the processes taking place in this regime 
is essential to obtain a complete picture of the process of 
accretion and the driving of outflows from low mass proto- 



stars as well as the early chemical and physical evolution 
of circumstellar disks. 

One of the most effective probes of warm dense gas is 
through emission in molecular ro-vibrational transitions. 
The most common bands readily available from ground- 
based facilities are the CO overtones and H2 fundamental 



bands around 2.2 fim (e.g. Reipurth & Aspin 1997) and 
the fundamental transitions of CO in the M-band around 
4.7 /im. Also emission from very hot water gas (~ 2000 K) 
near 2.29 /im has been reported toward a few sources (e.g. 



Najita et al. 200C). Since the upper levels of these transi- 
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tions lie at temperatures of up to a few thousand Kelvin, 
they probe hot gas with temperatures between 100 and 
1000 K, making them ideal to study the region of interac- 
tion between disk, protostar, accretion and outflow. With 
the new generation of sensitive ground-based high resolu- 
tion spectrometers for the 3 — 5 /im region an efflcient 
window has been opened for the detailed study of the 
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CO ro-vibrational lines toward embedded young stellar 
objects (YSO). The past generation of instruments was 
suitable to either observe a low resolution spectrum with 



a fairly wide spectral range (e.g. Teixeira et al. 1998) or a 
high resolution echelle spectrum with a very narrow spec- 
tral range (e.g. Carr et al. 2001). The main exception are 
the Fourier Transform Spectroscopy (FTS) observations 
of Mitchell et al. (1985, 1990), who observed a number 



of high-mass stars in the entire M band at high spectral 
resolution (i? > 10^), but such studies are limited to the 
brightest objects. VLT-ISAAC has a large instantaneous 
spectral range in the M-band (0.237 fim), which combined 
with a medium resolution of i? = 5 000 — 10 000 and a lim- 
iting magnitude of M ~ 10 allows the entire fundamental 
band of gaseous CO of a wide range of young stars to be 
observed in a short time, including low mass stars down 
to a few tenths of a solar mass in the nearest star-forming 
clouds. 

The embedded stars studied so far in CO ro-vi brational 
bands have showed mostly lines in absorption ( Mitchelj 
et al. j990t [Boogert et al. 2000| , |2002| ) implying that the 
warm CO gas is seen in front of a bright infrared contin- 
uum, produced by hot dust close to the central object. The 
fundamental CO lines are usually only seen in emission 
towards more evolved sources characterized by a class II 
type spectrum where a circumstellar disk is directly vis- 
[Blake et al. 2002|). CO overtone 



ible (Carr ct al. 2001 



bandhead emission at 2 /im has been observed in emis- 
sion toward a few intermediate mass pre-main sequence 
stars ( [Thompson 1989| ; [Najita et al. 1996| ) and T Tauri 
stars (Garr 1989), and has been associated with hot gas 
(T - 1 500 - 5 000 K) located within a fraction of an AU 
in a Keplerian disk. 

We present here the peculiar 4.5 — 4.8 /Ltm spectrum 
of the illuminating source IRSl=Elias 21 of the bipolar 
reflection nebula GSS 30 located in the core of the p 
Ophiuchus molecular cloud at a distance of 160 pc. It has 
been classified as a low mass class I YSO from its spectral 
ener gy distribution (SEP ) ( |Grasdalen ct al. 1973| ; Elias 
1978 |; linking et al. 1989) and low bolomctric luminosity 
(Lboi ^21 — 26 L(^, [Greene et al. 1994 : Bontcmps et al 
200l[)~[Extensive polarimetric studies in the H and K band 
of the reflection nebula have shown that the source is sur- 
rounded by a large (~ 2 000 AU) disk-like envelope and 
a smaller circumstellar disk of ~ 150 AU, which are in- 
clined about 25° away from the plane of the sky, i.e. close 
to edge-on (Chrysostomou et al. 1996, 1997). The high 



degree of linear polarization (up to 50%) as well as the 
presence of circular polarization imply that the light com- 
ing from the reflection nebula must have been multiply 
scattered, before heading into the line of sight. 

Two other sources (IRS2 and IRS3) are present toward 
the JsT-band reflection nebula (within 20"of IRSl). IRS2 
has a class III SED and is probably a more evolved star 



IRS3 shows strong 6 cm emission (see Leous et al. 1991 
where IRS3 is designated LFAM 1). 

The presence of a molecular outflow has not been 
firmly established. High velocity red- and blue-shifted CO 
millimeter emission to the south of GSS 30 has previously 
been reported by Tamura et al. (199C| ), but since both 
lobes are located to the south of the infrared source, the 
gas is likely to be associated with the VLA 1623 jet, which 
is passing only 30" to the SW of IRSl. Using millimeter 
interferometric line data, Zhang et al. (1997) find evidence 
for a spherical expansion of the core surrounding the three 
sources in GSS 30. In addition, the inner region of the core 
seems to be cleared, which is indicative of a young outflow. 
This is additionally supported by the presence of variable 
unbound water maser emission from within 0'.'3 of IRSl 



( Claussen ct al. 1996 ) 



( [Andre fc Montmerle 1994[ ). IRS3 has a class I SED, but 
is much fainter than IRSl in the near-infrared. It has a 
bolometric luminosity of 0.13 Lq, (Bontcmps et al. 2001). 



2. Observations 

The observations of GSS 30 IRSl were carried out using 
the long wavelength (LW) medium resolution mode on 
the Infrared Spectrometer And Array Camera (ISAAC) 
mounted on the Very Large Telescope (VLT-UTl) at the 
Paranal Observatory on September 3, 2001. The detector 
for 3 — 5 /im observations is a IK x IK Aladdin InSb array, 
which allows a spectral coverage of 0.237 /um per setting 
in the Af-band with a spectral resolution of up to i? = 
10 000. The detector resolution in the spatial direction in 
the spectroscopic modes is O'.'148/pixel. The spectrum of 
this source was taken as a part of a large program to study 
ices and gas around young low mass stars (van Dishoeck 
[et al. 2001[ ). 

The data were obtained under excellent photometric 
conditions in service mode with 0'.'3 infrared seeing and 
< 10% humidity which gave a very stable atmosphere. The 
frames were both chopped and nodded along the slit with 
a 15" chop throw. The 0"6 slit was used in two settings for 
20 minutes per setting yielding a final signal to noise ratio 
on the continuum of the extracted spectrum ranging from 
50 to more than 100 over a spectral range from 4.5 /im to 
4.82 /im and with a spectral resolution of i? = 5 000. 

The standard star BS6084 (BlIII) was observed at air 
masses between 1.2 and 1.5 immediately before GSS 30 for 
removal of telluric features, while the source was observed 
at air masses between 1.5 and 2.1. Because the source 
was already descending at the time of observation, it was 
not possible to observe source and standard at the same 
airmass as is otherwise of importance in order to obtain a 
good telluric subtraction. However, due to the exceptional 
stability of the atmosphere through most of the night we 
were able to correct for the airmass difference using a sim- 
ple Lambert-Beer law, which states that the depth of the 
atmospheric lines scales exponentially with the airmass. 
The final corrected spectrum turned out to be of accept- 
able quality in spite of the airmass difference. Our expe- 
rience is that it is crucial to use a B star or an early A 
star as the standard because many narrow intrinsic stellar 
lines become visible in the 3 — 5 /xm region for spectral 
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Fig. 1. The VLT-ISAAC Af-band spectrum of GSS 30 IRSl at R = 5 000. The statistical (3ct) error is shown below 
the spectrum. Note that this error does not take the systematical errors introduced by variable telluric absorption 
into account which dominate in some parts of the spectrum. Since the peaks in the statistical error indicate deep 
atmospheric absorption they can be used to locate where telluric residuals might be interfering with the intrinsic 
spectrum. The telluric residuals which are exceeding the noise are marked with a © and are seen to correlate well with 
the peaks in statistical error. 



resolutions better than ^ 1000 if an F or later type star 
is used. 

The spectrum was reduced with IDL routines using 
standard methods appropriate for chopped and nodded in- 
frared spectroscopy. Bad pixels were removed, the frames 
were distortion corrected using the telluric features in a 
star-trace map as reference, and AB nodded pairs were 
subtracted to remove the long time scale differences which 
the chopping is not capable of removing. The spectra were 
extracted by integrating over a spatial width determined 
by the positions along the slit where the signal drops be- 
low the 3(7 level. The raw spectrum was then ratioed by 
the standard after applying an optimized shift and air- 
mass correction. In the final spectrum, the points where 
telluric residuals are apparent have been removed, includ- 
ing the region short-wards of 4.55 /im, where atmospheric 
CO2 lines render the spectrum useless. The final result is 
shown in Fig. |^. 

The CO ro-vibrational lines in the spectrum are red- 
shifted with ^ 20 km s^^ due to the systemic velocity of 
the Earth relative to the source, which means that there 
is a significant overlap between intrinsic CO lines and tel- 
luric absorption by CO. In spite of this, the telluric di- 
vision is very good, especially when taking the unfavor- 
able airmasses of source and standard into account. A few 



residuals from telluric ^ CO absorption are visible and are 
marked in the final spectrum. Since the resolution of our 
spectra is 60 km s^^, none of the features in the source 
are unaffected by the telluric lines. At higher resolution 
{R > 20 000) the telluric lines may saturate, thereby forc- 
ing a division by zero when ratioing with the standard 
and resulting in a loss of information in the affected pix- 
els. At sufficiently high resolution, this will not be a prob- 
lem since the narrow telluric lines will often be completely 
offset from the intrinsic lines. However, at a resolution of 
5 000, the typical telluric line has a depth which is only a 
fraction of the intrinsic continuum allowing a significant 
signal to be present in the standard spectrum even in the 
center of a telluric line and only very few lines are partially 
or completely lost due to saturated telluric lines. 

By measuring the strength of the small residuals from 
telluric lines of species not intrinsically present in the 
source spectrum, the systematic uncertainty in the source 
lines is estimated to be less than 15%. The statistical un- 
certainty for the brighest lines is typically a few %, while 
the faintest lines are 3cr detections. The spectrum was 
wavelength calibrated relative to the telluric lines and flux 
calibrated relative to the standard star. The wavelength 
calibration of high resolution M band spectra is very ac- 
curate due to the large number of telluric lines, and a 
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conservative estimate for the precision is ~ 5 km s^^ over 
the entire spectrum. Indeed, the standard deviation of the 
Hne center velocities given in Table|l|is4.2kms~i. Due to 
the difference in airmass between standard and source, we 
estimate the systematic uncertainty in absolute flux cali- 
bration to be ~ 30%. Since the bright lines are typically 
most severely affected by telluric residual, all line fluxes 
have total uncertainties with values between 15 and 50%. 
The exact uncertainty for each line is difficult to estimate 
and we consequently adopt a value of 30% for all lines. 

3. Results 

3.1. Signatures of hot CO gas 

It is evident that the A/-band spectrum of GSS 30 is dom- 
inated by strong emission lines from the ro-vibrational 
(w=l-0) transitions of ^^CO gas. R- and P-branch lines 
within the observed spectral range are detected up to 
J=15 and 7=18, respectively. The lines are all unresolved 
at a resolution of i? = 5 000, which gives an upper limit to 
the intrinsic line width of 30 km s^^. No apparent decrease 
in line intensity is seen towards higher rotational quantum 
numbers, so observing the source in the entire A/-band 
will clearly detect more lines. The brightest lines reach 
peak fluxes of more than 200% of that of the continuum. 
However, since none of the lines are resolved, this is only a 
lower limit to the peak flux. Fainter detected features in- 
clude i^CO v=l-0 lines as wefl as ^^CO v=2-l lines. The 
fact that the observed gas has a population excited to the 
second vibrational level already indicates that the gas is 
at least warmer than a few hundred K. The ^■^CO ■i;=l-0 
lines have intensities of order 2 x 10^"'^'^ erg s^^ cm~^, the 
13C0 u=l-0 and ^^CO v=2-l lines being 1-2 orders of 
magnitude fainter. The line fluxes of all unblended lines 
were extracted by fitting a gaussian as well as a first order 
polynomial to the local continuum. After subtraction of 
the continuum, the fiuxes were derived by integrating the 
spectrum in a 3cr region around the line. 

The extracted fluxes of all unblended lines are given 
in Table [|. Many of the faint lines are completely blended 
with the bright lines, but it was still possible to extract 
12 I^CO v=l-0 fines and 4 ^^CO ■i;=2-l lines, which is 
sufficient to constrain the most important parameters of 
the emitting gas. 

In addition to the CO emission lines, a faint 4.67 fiia 
CO ice feature is visible. Since the feature is partly filled 
by the bright ^^CO P(l) fine it is not possible to obtain 
a well-defined ice profile. For the same reason the affected 
emission line is not included in the analysis and model- 
ing. There is also a very broad (FWHM ^ 500 km s^^), 
4.65 fiia Pfund /3 hydrogen recombination line visible, al- 
though with the bright CO lines superposed. 

Finally the II2 — S(9) line is detected with an inte- 
grated line fiux of 2 ± 1 x 10""'^'* erg s~^ cm~^, signifying 
the presence of hot (> 1000 K) molecular hydrogen or 
possibly shock-excited gas. 



Table 1. Integrated CO line fluxes toward GSS 30 IRSl 
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This line is severely affected by telluric residual. 



3.2. Size of the emitting region 

In order to determine the parameters of the line emitting 
gas unambiguously, the size of the emitting region must 
be known. The FWHM of the continuum in the spatial 
direction is about 0'.'5, which is consistent with the seeing 
throughout the observation. However, the lines appear to 
be spatially resolved, although most of the line flux is cen- 
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Fig. 2. Spatial cross section of the resolved line emission. 
The spatial distribution of the emission (along the slit) is 
shown averaged over all ^^CO v=l-0 lines with the con- 
tinuum spatial profile subtracted. The line spatial profile 
was normalized to the continuum profile to show that the 
lines are detected over a larger region than the continuum. 
The direction to the north on the sky is indicated on the 
figure. 



tered on the source. The wings of the spatial profile of each 
line extend up to 2" away from the central source along 
the slit. This is illustrated in Fig. |2[ which shows the nor- 
malized average of the spatial profile of the lines, with the 
normalized average of the continuum profile subtracted. It 
is evident that the line emission is resolved and detected 
to a distance of more than 320 AU from the central source. 
There is also an indication that the northern part of line 
emission is stronger than the southern part, which agrees 
with the morphology seen in the L-band image of IRSl 
shown in Fig. |^. 

This extension suggests that the line emission is associ- 
ated with the scattering nebula and is not coming directly 
from the central source. Since the continuum is less spa- 
tially extended than the line emission, the line emitting 
region appears to be spatially distinct from the continuum 
emitting region. Possible interpretations are that the line 
emission can be coming from gas present in the bipolar 
lobes (e.g. an outfiow), from resonance scattering of con- 
tinuum emission or from scattering of lines emitted closer 
to the star. Each of these possibilities will be discussed 
in detail below as the size of the emitting region depends 
strongly on the scenario. 

4. Models 

4.1. Optically thin LTE models 

As a simple first approximation, the extracted line intensi- 
ties were fitted with a one-temperature optically thin LTE 
model. The molecular data were taken from the HITRAN 
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database (Rothman et al. 1998) and the relative line in- 
tensities were calculated directly from a Boltzmann level 



Fig. 3. The 3.21 fim acquisition image of the reflection 
nebula surrounding GSS 30 IRSl. The north-eastern lobe 
is clearly visible and also extended emission from the 
south-western lobe is resolved. The slit was centered on 
(0, 0) and directed north-south. The contours are roughly 
logarithmic with the lowest contour at values two orders 
of magnitude smaller than the highest contour. 

population distribution. The two parameters of the model, 
the gas kinetic temperature and the CO column density, 
should be well constrained in the optically thin limit. If 
only the ^^CO 'y=l-0 transitions are used, an optically 
thin model is to some extent degenerate in the two pa- 
rameters. This degeneracy is due to the fact that the pop- 
ulations in the vibrationally excited levels are much more 
sensitive to temperature changes in the 400 — 800 K range 
than the rotational level populations. Since lines from two 
species are available, the degeneracy is broken. It is rea- 
sonable to assume that the emission from both species is 
coming from the same mass of gas, and tha t the ^^C/^'^C- 
ratio is close to 60 (e.g. Bensch et al. 2001 ). 

This model is not able to provide a good fit to the 
entire line spectrum. However, from Fig. ^ it is evident 
that the faint ^^CO w=l~0 and ^^CO v=2-l lines can be 
simultaneously fitted with a single temperature gas in the 
optically thin limit. The model lines of ^^CO v~l-Q are 
however an order of magnitude too bright. Excluding the 
^^CO t;=l-0 lines from the fit yields a well-constrained 
temperature of 515 ± 25 K and a total ^^CO column den- 
sity of 4 ± 1 X 10^** cm~^. Since the fit includes lines from 
both CO isotopic species, we conclude that it is unlikely 
that the apparent discrepancy between the ^^CO 'i;=l-0 
and ^"^CO v—l-Q lines is due to a deviation from the nor- 
mal interstellar ratio. The simplest explanation is that the 
^^CO v—l-Q lines are optically thick, while the rest remain 
optically thin. This would conveniently prevent us from 



K. M. Pontoppidan et al.: Bright CO ro-vibrational lines in GSS 30 IRSl 




4.58 



4.60 4.62 
Wavelength [/xm] 



4.64 



1.0 
0.8 



^ R13 



0.6 — 



0,4 
0.2 



K12 



R11 



0.0 

-0.2 b 




4.66 



4.68 



4.70 

Wavelength [fim] 



4.72 



4.74 



1.0 
0.8 
" 0.6 



X 



0.4 
0.2 
0.0 
-0.2 



1-|_ P|4 P|5 



p.i 



P.6 



P,2 



P7 



T 

P.8 





4.76 



4.78 



4.82 



Wavelength [/xm] 



4.84 



Fig. 4. Continuum subtracted spectrum of GSS 30 IRSl. Top row of line identifications show the ^'^CO ■i;=l-0 lines. 
Bottom row of identifications show the ^^CO v=2-l lines. The data have been over-plotted with the optically thin 
LTE model spectrum. The model ^^CO v—1-0 line fluxes have been divided by a factor of 7 for clarity and to illustrate 
the optical depth effects. 



seeing the whole column of ^^CO and could possibly lower 
the line intensities enough to explain the data. However, 
the fundamental ro-vibrational transitions of ^^CO do not 
produce optically thick lines for column densities less than 
^ 10^^ cm~^, which is still 250 times more than the sim- 
ple model dictates. Thus the optically thin LTE model is 
internally inconsistent. 

4.2. Results of the full radiative transfer model 

In order to quantify the discrepancy more accurately, we 
computed the full radiative transfer for a single temper- 
ature static envelope in LTE. For simplicity a spherical 
geometry was assumed since the results are not very sen- 
sitive to geometry. The size and density of the emitting 
region cannot be constrained by the model, so the quanti- 
ties fitted were the single gas temperature and the column 
density of the gas, by keeping the size fixed at the observed 
4". 

Furthermore, it is assumed that the line intensity is 
independent of the underlying continuum. This last as- 
sumption may not be completely valid, since the CO lines 
could be present in absorption in the continuum before the 



lines from the emitting gas are superposed. The effect of 
this would be to lower the line peak flux with maximally 
the continuum flux level. Since the lines are completely 
unresolved, the true peak flux of the lines is higher than 
that observed. The observed peak flux of the ^^CO v=l-0 
lines is 2-3 times higher than the continuum flux, so we 
conclude that this effect can only effectively lower the in- 
tensity of these lines with a small fraction at most. The 
fainter lines can be significantly suppressed by absorption 
in the continuum if they are optically thick, but this would 
only increase the discrepancy between the bright and faint 
lines. Also, if the absorption lines are significantly broader 
than the superposed emission lines, the observed line in- 
tensities could be sufficiently suppressed. However, since 
we have no possibility of constraining this with the present 
data, we assume that the possible absorption lines and the 
emission lines have similar intrinsic widths, i.e. the con- 
tribution from absorption in the continuum is negligible. 
Higher resolution spectra may clarify the matter. 

The ^"^CO v—l~0 lines were first fitted independently 
of the i^CO and ^^CO v=2-l lines. The model 

confirms that the two sets of lines cannot be fitted si- 
multaneously, although each set is consistent with a sin- 
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Fig. 5. maps of the best fitting models. Left paneh the ^^CO v—1-0 and ^^CO v—2-1 hnes plotted against the 
^•^CO v=l-0 lines without the assumption of scattering, in which case the sets of well-fitting parameters do not overlap. 
Middle panel: contours for the ^^CO v—1-0 and ^^CO v=2-l lines over-plotted by contours for the ^'^CO v=l-0 
lines, but now with a scattering factor e = 1.4 x 10'^. Right panel: The weighted sum of the two maps from the 
middle panel. The contours show in all the maps the 1, 2 and 3<7 levels. 



gle temperature gas. There is a certain degeneracy in the 
two free parameters, especially in the ^'^CO lines, and 
the model gives a one-dimensional family of good fits for 
each species. However, the two families do not overlap 
for the observed line intensities. Fitting the ^^CO ■i;=l-0 
and the ^^CO v=2-l lines simultaneously yields poorer 
fits than fitting them separately. All the lines are in 
the extreme optically thin limit with typical values of 
T ~ 10"'^ for the "'^^CO v=l-0 lines. The discrepancy be- 
tween the ^^CO t;=l-0 and "'^^CO v=2-l lines may be due 
to non-LTE effects, reflecting that the rotational temper- 
ature Trot and the vibrational temperature Tvib are not 
equal. However, non-LTE effects cannot explain why the 
I^CO v=l-0 to i^CO line ratios differ with almost 

an order of magnitude from the predictions of the simplest 
LTE model, since the excitation structures of different iso- 
topic species must be almost identical. 

To explain the data, we propose that all the observed 
line emission has been scattered by dust in the bipolar 
cavity seen in the near-infrared. The scattering efficiency 
is very low at 4.7 fim, and thus only a small fraction of 
the original line emission will be scattered into the line 
of sight. This means that the line intensities from the 
emitting region and consequently the column densities can 
reach much larger values than those previously suggested 
by the model. To test this hypothesis, we adopted a factor, 
e = 5* X {TrR'^)/n, which accounts for our lack of informa- 
tion about the true solid angle of the emitting region Q, 
and the ratio of intrinsic line flux to line flux scattered 
into the line of sight, S. R is the observed angular ra- 
dius of the emitting region, which is taken to be i? = 2". 



with this factor. Multiplying all observed line intensities 
with a sufficiently large factor will allow the ^^CO v=l-0 
lines to become optically thick. We are able to constrain 
the scattering factor in the model by demanding a simul- 
taneous fit to all the lines. 

The best fitting value of the scattering factor is e = 
1.4 X 10'^ with a reduced of 1.29, which gives a well- 
constrained single temperature LTE gas at 515 ± 5 K and 
N(CO) = 1.5 ± 0.5 X 10^^ cm-2, as indicated in Table § 
With an extent of the emitting region of 4" an upper limit 
to the total gas mass of 200 M® is found. Since the total 
gas mass scales with $7°-^ an absolute lower limit to the 
mass is 1.3 M0 with a diameter of the emitting region of 
0.19" = 30 AU, assuming that S* = 1 and ignoring that 
line emission is observed to much larger distances from 
the central source. Since S must be significantly less than 
unity, the total mass is most likely significantly larger than 
the lower limit. It is not necessary to invoke any non-LTE 
effects to obtain a nearly perfect fit to all the observed lines 
and the assumption of LTE holds such that Tjot = ^vib. 
The maps of the best-fitting models are shown in Fig. 



It is possible to make an estimate for an absolute lower 
limit to the grain size by assuming = 1. In this case the 
scattering efficiency, Qsca > since it is not possible to 
scatter more light into the line of sight than out of it. This 
happens for silicate grain sizes larger than 0.2 — 0.3 fj,m 
(Draine & Lee 1984; Draine 1985). As a consequence the 



grains must be large, perhaps as large as 1 /im, which is 



expected for grain growth in dense media (Pendleton et al 



Absolute upper and lower bounds to the size of the emit- 1990). On the other hand, the lower limit on the size of the 



ting region are given by the natural conditions S > 1 and 
< ttR^. Decreasing the size of the emitting region will 
increase the column densities, but since line emission is 
observed far from the source, we cannot explain the entire 
discrepancy with a smaller emitting size. To recover the 
original ffux it is necessary to multiply the observed values 



scattering grains docs not fit with the polarization maps of 



the GSS 30 nebula. Chrysostomou et al. (1996) finds that 
grains larger than 0.35 /im are inconsistent with the vari- 
ation of linear polarization with wavelength. However, the 
polarization was mapped on a much larger angular scale 
than the size of the nebula at 4.6 /zm, so the discrepancy 
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Table 2. Best fitting models for GSS 30 IRSl 



Model 



Lines used 



^^CO v=l-0 + '"CO v=2-l 
''■^CO v=l-Q 

^^CO v=l~0 + ^^CO u=2-l 
All 



T 
[K] 



N(CO) 



Red. 



Optically thin direct calculation 

Full radiative transfer 

Full radiative transfer 

Full radiative transfer with scattering 



500 ± 25 
degenerate 
725 ± 50 
515 ±10 



4±lxl0" 1.4 

degenerate 0.85 

2.5±lxl0" 2.38 

1.5 ± 0.5 X 10^** 1.29 



may be explained by grain growth within a few hundred 
AU of the central source. 



4.3. Other scenarios and pumping mechanisms 

One way of avoiding the scattering assumption entirely 
may be to apply a two-temperature model to the data. In 
this scenario a small part of the gas close to the central 
source with a kinetic temperature of ~ 1000 K would be 
responsible for the ^^CO v=2-l lines as well as part of the 
high J lines of the i;=l-0 lines, while a massive outer com- 
ponent at ~ 200 K would dominate the emission in the low 
J lines. The scattering assumption would not be needed 
since the ^^CO v=l-0 lines require a much higher column 
density to create the observed line intensity at 200 K com- 
pared to 500 K and will thus more readily become optically 
thick. We applied a two-temperature model by requiring 
that the total column density must be the same as that 
determined by the single-temperature model. This is rea- 
sonable because the optical depth of the lines is given by 
the ratio between the ^^CO and ^^CO lines independently 
of e and therefore must be conserved. As in the single- 
temperature model three parameters are varied, namely 
the two temperatures and the ratio of the column density 
between the components. In other words, e is exchanged as 
a free parameter with a second temperature, thus keeping 
the number of free parameters constant. 

The "'^^CO lines can be fitted well within the two- 
temperature model, although the parameters are some- 
what degenerate in the same way as in the single temper- 
ature model. However, it was not possible to fit the ^^CO 
lines within this scenario, with the best reduced being 
of order 10. Furthermore there is a preference towards a 
single low temperature in the ^^CO lines which is incon- 
sistent with the ^^CO lines. Inclusion of all lines gives no 
meaningful fits within the parameter space. The explana- 
tion is that it is not possible to combine the ^^CO v—2-1 
and the ^^CO v=l-0 lines, since the high temperature 
component easily produces too much ?;=l-0 emission in 
order to fit the t;=2-l emission. We conclude that a two- 
temperature model is clearly inferior to the scattering 
model in explaining the observed line emission. 

Radiative pumping by a strong infrared continuum is 
another possibility that needs to be considered as an al- 
ternative excitation mechanism and a way to avoid the 
scattering scenario by introducing non-LTE effects. This 



would produce a single vibrational temperature (3coville 
et al. 1980). However it does not explain why the rota- 



tional and vibrational temperatures seem to be identical. 
Also, luminosities larger than 10'^ L0 are required to ef- 
fectively pump CO molecules to the second vibrational 
level or higher at distances from the central source be- 
yond 1 AU ( ^coville et al. 1980| ). Using the 1-dimensional 
Monte-Carlo radiative transfer code from Schoier (2000[ ) 
illustrates the effects of radiative excitation on the line 
populations. The simplest possible model for radiative ex- 
citation assumes a blackbody with a luminosity L — 25 Lq 
at the center of a single density spherical and static enve- 
lope . The blackbody temperature is taken to be the colour 
temperature of the continuum around 4.7 /im, determined 
from the M-band spectrum to 420 K. The radius used is 
the observed 320 AU, the gas kinetic temperature is as- 
sumed to be a power law beginning at the surface of the 
central source. The molecular constants are taken from 
^choier et al. (2002 ). This leaves the gas density as the 
only free parameter, which basically determines the opti- 
cal depth of the lines and consequently the importance of 
line trapping effects. The model calculations show that a 
density of at least riHa ~ 10^ cm~'^ is needed to reproduce 
the approximate flux level of the ^^CO w=l-0 lines. The 
vibrational excitation temperature is largely determined 
by the continuum flux level at 4.7 fim and therefore by 
the luminosity and to a smaller extent the temperature 
of the radiating source. Since the dominating mechanism 
for populating the v=2 level is via the v=l level, and 
since both the w=2-l and 11=1-0 transitions are placed in 
the same narrow wavelength region, the radiative excita- 
tion rates are strongly dependent on local photon density 
in the 4.7 /zm region. Consequently, as long as the di- 
rect route through the v=2-0 overtone transitions around 
2.35 fxm is not important, the vibrational level population 
depends on the temperature of the input radiation field 
only through the flux level at 4.7 /im. At the adopted lu- 
minosity, the ^^CO v=2-l line fluxes are two orders of 
magnitude too low, reflecting an excitation temperature 
of 400 K. To significantly populate the second vibrational 
levels, a luminosity of at least 10^ is required, in agree- 
ment with Scoville et al. (1980). Finally, the population of 
the rotational levels shows significant departures from a 
Boltzmann distribution, especially when line trapping be- 
comes important. Altogether, this indicates that radiative 
pumping as an excitation mechanism in the case of GSS 
30 IRSl is unlikely. 

Finally, the possibility that the lines are continuum 
emission which is resonance scattered in CO transitions 
should be considered. It is known that this takes place in 
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the case of asymptotic giant branch stars, e.g. in the case 
of Mira ( |Ryde fc Schoier 200lD . However, since this mech- 
anism probes the scattering gas rather than the emitting 
medium, it will not solve the discrepancy. Also it would 
require that the gas at a distance of 320 AU from the cen- 
tral source has a single temperature of 500 K, which seems 
highly unlikely. 



4.4. Hydrogen recombination lines and mass loss rate 

There are a number of hydrogen recombination lines com- 
mon to YSOs (Pf/3, Bra, Pf(5) present in the 3 — 5 fim re- 
gion, which are usually taken as an indicator of accretion 
activity. This is particularly true in the case of low mass 
young stars, where the hard radiation from the boundary 
layer between disk and star is needed to produce the ion- 
izing radiation, since the central star is not hot enough. 
Because these lines have broad wings (3 — 500 km s^^), 
they are often interpreted as coming from a strong stellar 
wind, which may be the engine of an outflow. The possi- 
bility that the CO gas could be associated with a wind or 
molecular outflow should therefore be explored. If associ- 
ated with a stellar wind, the lines could be formed either 
as cooling emission from the gas as molecules are reformed 
or as emission from a shock as the wind collides with the 
circumstellar envelope or disk. It is known that no large 
scale outflow is present which could be produced by any of 
the three stars in the cluster, but there is some evidence 
of an expanding core surrounding the cluster from inter- 
fcrometric observations of rotational ^'^CO and C-'^^O lines 



(Zhang et al. 1997) 



Hydrogen recombination lines emitted from a plasma 
moving at supersonic speeds can be modeled in the 
Sobolev approximation (Sobolev 1960), which is valid 
when the thermal line broadening is much smaller than the 
wind speed. This formalism has been treate d for spherica l 
ioni zed winds in LTE from evolved stars in Castor (1970 ) 
and Krolik fc Smith (1981 ), and for Herbig AeBe stars in 
Nisini et al. (1995 ). Using this formalism an ionized mass 



loss rate can be estimated for the star if a wind origin for 
the hydrogen lines is assumed. The extinction corrected 
line intensities for the hydrogen recombination lines ob- 
served in GSS 30 IRSl are tabulated in Table [s] along 
with the mass loss rates derived from each line. The line 
ratios and the mass loss rate are not sensitive to the tem- 
perature, so a constant electron temperature of 10* K has 
been assumed. Finally, it is assumed that the central star 
does not shadow a signiflcant part of the radiation from 
the wind. Changing the model parameters, i.e. the tem- 
perature and the size of the ionized region, with a factor 
of two changes the mass loss rate with a factor of two at 
most. The velocity structure is such that the wind starts 
at the surface of the star with 20 km s^^ and rises quickly 
to 450 km s^i to fit the observed FWHM of the Pf/3 fine. 
The exact shape of the velocity structure is otherwise not 
significant for the derived mass loss rate. 



We find that an ionized mass loss rate of 1.7 x 
10^^ Mq year^^ = 0.6 M0 year^^ is consistent with the 
two observed Pfund lines. The Bra line is too faint for this 
mass loss rate, which may indicate deviations from LTE. 
We can rule out that it is an effect from a poor extinction 
correction since the Bra line lies between the Pf/3 and PiS 
lines in wavelength. The most serious problem with this 
model is likely that very little is known about the structure 
of winds from low mass embedded YSOs; in particular, 
the ionization structure will depend directly on energetic 
processes which are poorly understood, such as accretion 
activity or interactions of the circumstellar matter with 
strong magnetic fields. Also, if the wind is predominantly 
neutral, then the total mass loss rate will be much higher 
than the ionized mass loss rate and will become inconsis- 
tent with the non-detection of an outflow. It is possible 
that other mechanisms such as accretion flows must be 
considered for the formation of recombination lines from 
low mass young stars. The hydrogen lines are therefore 
not likely to provide any strong constraints on the emit- 
ting CO gas and a detailed discussion on the hydrogen 
emission is outside the scope of this paper. 

5. Discussion 

In summary, to determine the origin of the CO gas emis- 
sion it must be explained why only a single well-defined 
temperature is needed. Which mechanism can heat up to 
100 of gas to a unique temperature of 500 K, yet 
keep the intrinsic line width less than 30 km s^^? Finally, 
why is no other similar embedded source from the litera- 
ture showing the same strong CO rovibrational emission 
as GSS 30 IRSl? 



5.1. Outflow? 

If the warm CO gas resides in the wind/outflow compo- 
nent of the circumstellar environment then an estimate of 
the time needed to create it can be found from the mass 



loss rate derived in Section 4.4 under the assumption that 



the wind is predominantly ionized. A mass between 3 and 
100 A/0 corresponds to a production time of minimally 
5 years and maximally a few centuries, while the most 
likely time is about one century. It seems unlikely that 
hot, ionized material emitted over a longer period should 
thermalize at a single temperature. Also, an outflow ori- 
gin of the line emission would produce broad wings in the 
lines, which is not observed. If the wind is predominantly 
neutral, the resulting mass loss rate would be so high that 
a clear outflow should have been detected. 



5.2. Inner disk? 

Another possibility is that the emission is produced by 
warm thermalized gas in the disk itself. Since more evolved 
T Tauri disks are known to exhibit similar CO emission, 
although with smaller intensities, it is conceivable that we 
are seeing an equivalent process in the case of a younger 
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Table 3. Hydrogen recombination lines observed in GSS 30 IRSl 



Transition 


Line intensity 




FWHM 


Mass loss rate 




[erg cm^^ s^^] 




[km s^^] 


[Mq year-i] 


Pf/3 


1.3 ±0.5 X 10" 


13 


450 ± 50 


1.7 X 10"" 


Bra 


1.6 ±0.2 X 10" 


13 


< 500 


0.8 X lO"*^ 


PiS 


6.0 ± 1.0 X 10" 


14 


< 500 


1.8 X 10"® 



disk. A typical circumstellar disk around a low mass young 
star is often observed to have a mass of a few times the 
minimum solar nebula, i.e. Mdisk ~ a few 0.01 Mq (Andre 



& Moi jtmcrle 1994| ; [Osterloh fc Beckwith 1995| ). The disk 
mass for GSS 30 IRSl inferred by the 1.3 mm continuum 
emission is Mdisk = 0.03 Mq ( [Andre fc Montmcrlc 1994| ). 
In hydrostatic equilibrium, a disk irradiated by the central 
star has a surface density E ~ R~^/^ (Chiang & Goldreich 
1997l3lf viscous dissipation in the disk is taken into ac- 
count the surface density attains a flatter i?-dependency, 
E ~ R-^. In both cases the accumulated disk mass reaches 
the observed 10 — 100 between 2 and 10 AU from the 
central star. A single temperature of 515 K is however not 
consistent with the disk models, which prescribe a large 
range in temperatures with a small 1000-2000 K compo- 
nent within a few stellar radii to an extended 100-200 K 
component within a few AU of the central star. However, 
in no models for passive circumstellar disks is it expected 
that the disk exhibits a single temperature, indeed quite 
the contrary is the case. Furthermore, Keplerian rotation 
within a few AU would produce observable line broaden- 
ing. The narrow lines could only be explained if we are 
seeing light emitted vertically from the plane of the disk 
before being scattered into the line of sight, effectively re- 
moving the effects of rotation. 



5.3. Accretion shock? 

A final option is that the emission lines are cooling lines 
from post-shocked dense gas. Vibrational H2 emission 
is one of the principal tracers of low density (wh- ^ 



10 cm shocked gas. As mentioned in Section 3.1 the 



(0-0) S(9) line of H2 at 4.695 /im is seen in our M-band 
spectrum, while no H2 lines are seen in the if-band spec- 
trum of GSS 30 IRSl by [Greene fc Lada (1996| ) affhough 
no upper limit is given. Assuming that the observed H2 
(0-0) S(9) line is thermalized at the CO temperature, the 
required molecular gas mass is 5 M®. This is assuming 
that the H2 emission is seen directly and is not corrected 
for extinction, since the extinction is hard to estimate 
for embedded stars. A typical Ophiuchus extinction of 
Ay = 25 mag (Teixeira fc Emerson 1999) will increase 



the molecular gas mass by a factor of 2. If the emission in 
the H2 line is scattered in the same way as the CO emis- 
sion all values must be multiplied by e. If the H2 emission 
is thermalized at 515 K and directly observed then the 
integrated line flux expected for the 2.12 ^m (1-0) S(l) 
line is 7.5 x 10"^'^ erg cm~^ s~^, which should be observ- 
able. If the H2 lines are scattered into the line of sight 



with an efficiency of e, then no H2 lines should be visible 
in the A/-band. Consequently, the H2 line does not seem 
to be directly associated with the same gas emitting the 
CO lines since it is at least a factor of 10 too bright, but it 
may be an indication of a warmer component or the line 
may be pumped. Sensitive observations of H2 lines in the 
K and L band are needed to unambigously determine the 
origin of the molecular hydrogen emission. 

I ~l 

Neufeld & HoUenbach (1994) show that if a dense gas 

of 10^-5 < riHs < 10^2 cm is shocked, the vibrational 
emission from H2 is significantly suppressed partly be- 
cause the shock dissociates the H2 molecules and partly 
because the reforming molecules can easily be collisionally 
de-excited in the dense post-shock gas. The main coolants 
will be rotational and vibrational lines from H2O, OH 
and CO. This is valid for a wide range of shock velocities 
b <Vs < 100 km s^^ as long as the shock is of jump-type. 
A J-shock occurs in any case for > 30 km s"'^ and for 
smaller shock velocities if the magnetic field is weak or 
if the length scale over which the shock occurs is small 
compared to the length scale for acceleration of charged 
particles. Furthermore, the theoretical post-shock temper- 
ature structure gives an interesting prediction. A J-shock 
will dissociate the gas at the shock front. The gas then 
cools slowly from 10^ K principally through the Lya-line 
until CO forms at around 7000 K and iV(HI) - lO^i cm^^. 
Once the CO molecules are available the gas cools quickly 
to about 500 K through the CO vibrational lines, where 
a constant temperature plateau is maintained due to the 
release of chemical potential energy as H2 is formed which 
is balanced against the thermalized vibrational emission 
from molecules. The plateau occurs at an atomic hydrogen 
column density of 10^^ — 10^^ cm~^, which in the model 



by Neufeld fc HoUenbach (1994 ) corresponds to a CO col- 
umn density of 3 x 10^^ - 3 x 10^^ cm . It is stressed 
that the slower C-shock scenario was not treated in such 
high density models, and it may produce similar cooling 
lines. Further modelling is needed to exclude a slow shock. 
It is evident that the model values are remarkably similar 
to the values observed in GSS 30 IRSl with temperature 
Tgas = 515 K and column density iV(CO) = 2x 10^** cm^^. 
However, it is unclear why no hot 5000 — 7000 K com- 
ponent in the fundamental CO band is observed, which 
is also expected from the shock. Also no overtone emis- 
sion is seen in the i^-band spectrum of Greene fc Lada 
(19961 ). One plausible explanation is that the accretion 



activity is highly variable and that the hot component is 
only seen as the shock is occurring while the warm ther- 
malized gas is seen in the post-shock medium. This would 
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be supported by the variable water maser emission. The 
single observed H2 line is not inconsistent with an accre- 
tion shock although it is too bright to be associated with 
a dense shocked gas and may only be an indication of a 
warmer less dense component. 



As described in Ncufeld & HoUcnbach (1994), the most 
likely site for a dense shock in an embedded YSO is in a cir- 
cumstellar disk which is accreting matter. When infalling 
matter is colliding with the disk at supersonic velocities, a 
J-shock is expected to occur at the surface of the disk. If 
the accretion velocities are too low to produce a J-shock, 
the molecules will not be dissociated and the observed 
temperature plateau will not be present. The main un- 
certainty is if the high infall velocities required to create 
a dissociative shock are possible. However, infall veloci- 



ties of 10 km s within 10 AU are reasonable (Cassen 



& Moo rman 1981 ). Therefore, if the dissociative accretion 



shock scenario holds, the intrinsic line widths are predicted 
to be larger than or very close to 10 km s~^. An accretion 
shock at a large distance (> 10 AU) from the central star 
is consistent with the observation that the line emitting 
region is physically separated from the continuum emit- 
ting region. A shock will not heat a large amount of dust 
and the continuum emission, which has a colour tempera- 
ture of about 700 K, is likely to be produced by irradiated 
dust within 1 AU of the star. 

Other observable tracers of a shock in a dense medium 
are the rovibrational lines from the H2O 1^2 bending mode 
around 6 fim as well as the vibrational transitions of OH, 
all of which are unavailable from ground-based instru- 
ments. However, they may be detected by the upcom- 
ing Stratospheric Observatory For Infrared Astronomy 
(SOFIA) mission. 

The possible dissociation of the molecules has impor- 
tant consequences for the chemistry in the inner disk. If 
the emission is indeed coming from an accretion shock, 
and if a significant amount of matter passes through the 
shock, it will significantly affect the models of chemical 
evolution of the innermost parts of disks around low mass 
YSOs, by evaporating ice mantles and inducing high tem- 
perature chemistry. 



6. Summary 

In this paper, we have analyzed the 4.5 — 4.8 /im spectrum 
of the embedded young stellar object GSS 30 IRSl in the 
p Ophiuchus core. The observed emission lines from the 
rovibrational transitions of CO are fitted with a simple 
full radiative transfer ID model. It is found that the line 
emission must be scattered on a bipolar cavity in order to 
simultaneously account for the size of the observed emit- 
ting region, the absolute flux level of the lines and the 
ratio of ^^CO to '^'^CO lines. In this case the emission is 
well fitted by a single temperature gas with T=:515±5K 
and a column density of iV(CO) = 2 ± 0.5 x 10^* cm"^ 
Furthermore, assuming a two-temperature distribution of 
the gas does not yield satisfactory fits. 



The observed emission line spectrum can be best ex- 
plained by the presence of a J-shock in a dense medium, al- 
though other possible scenarios cannot be ruled out. Most 
likely the shock is produced by the accretion of gas onto 
a dense disk within a few tens of AU from the star as pre- 
dicted by Neufeld & HoUenbach (1994). Only this scenario 
explains why a large amount of CO gas is thermalized at a 
single temperature of ~ 500 K and why the emission lines 
are so narrow. 
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